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Abstract 



Based on the results in the previous papers that the boundary value 
problem y" — y 1 ' + y = y 3 , y(0) = , y (oo ) = 1 with the condition y (x) > 

for < x < oo has a unique solution y*(x), and a* — y* (0) satisfies 
< a* < 1/4, in this paper we show that y"~y'+y — y 3 , — oo < x < 0, 
with the initial conditions y(0) — 0, y'(0) = a* has a unique solution 
by using functional analysis method. So we get a globally well defined 
bounded function y*(x),— oo < x < +oo. The asymptotics of y*{x) as 
x — > — oo and as x — > +oo are obtained, and the connection formulas 
for the parameters in the asymptotics and the numerical simulations 
are also given. Then by the properties of y*(x), the solution to the 
boundary value problem r 2 f" + f = / 3 ,/(0) = 0, /(oo) = 1 is well 
described by the asymptotics and the connection formulas. 



1 Introduction 

In this paper, we use the results obtained in Q |J to study the boundary 
value problem 
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r 2 f" + f = f, 
f(r) - 0, as r - 0, 



< r < oo 



(1.1) 
(1.2) 
(1.3) 




If we make a transformation r = e x , f(r) = y(x), the equation is changed to 

y" -y' + y = y 3 - 

The results in 0] || can be briefly summarized as follows. The following 
boundary value problem 

II ' i 3 rv ^ ^ 

y -y+y = y,o<x< co, 

(P+) y(0) = 0,y(oo) = l, 

y(x) > 0, < x < oo, 

has a unique solution y*(x) for < a; < oo. The formula of 

a* = /(0), (1.4) 

is obtained, and 

< a* < 1/4. (1.5) 
In this paper, we show that the initial value problem 

y" -y' + y = y 3 , -oo < x < o, 

(P-) 2/(0) = 0, 
y'(0)=a*, 

has a unique solution by discussing the corresponding integral equation on 
a Banach space. Then we get a global solution y*{x) of y" — y' + y = y 3 for 
— cxd < x < +oo and the asymptotics as x — > —oo and as x — > +oo. And any 
bounded solution y(x) of this equation satisfying y(—oo) = 0, y(+co) = 1 is 
expressed as y(x) = y*{x — r), where r is the largest zero of y(x). Therefore 
we have solved the boundary value problem (|L^) ([L2; ) (|L^). 

This paper is organized as follows. In Sect. 2, y*(x) is extended to the 
negative axis by using contraction mapping theorem. In Sect. 3 we give 
numerical approximations for some important values including a*, which 
are used to represent asymptotics of the solution. In Sect. 4, the asymptotic 
expressions and connection formulas are given. 

2 Global Solution 

In ||], we have proved problem (P + ) has a unique solution y*(x) = y(x, a*). 
In this section we want to extend the solution to the negative axis. Consider 



the following problem 

y"-y' + y = y 3 , -oo < x < o, (2.1) 

(P-) i/(0) = 0, (2.2) 

y'(0)=a*. (2.3) 
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Let 

t = -x, u(t) = e~ x l 2 y(x). 
Then problem (P~) is equivalent to the following integral equation 
2 * . >/3, , 2 ft ^3 

where i > 0. Define 



u(t) = — ^=a*sin — t + -^= [ e' s sin — (t - s)u 3 (s) ds, (2.4) 



X = {f(t)\f is continuous on [0, oo), |/| < ~,t G [0,oo)}, (2.5) 

d(/,<?) = sup 1/(0-5(01, (2-6) 
te[o,oo) 

for f,g £ X. It is not difficult to prove the following lemma. 
Lemma 1 (X,d) is a complete metric space. 
Define 



s)u 3 (s)ds, (2.7) 



T(u)(t) = --j=a* sm + — ^ e " s sin ^-(t 
for u £ X. 

Theorem 1 T has precisely one fixed point u* in X. 

Proof. Let us first prove T(u) £ X when u £ X. By ( |2.7| ), T(u)(t) is 
continuous on [0, oo). By (|1.5|) , 

1 



which implies that if |u| < 1/2, 



So T{u) £ X. 

Next we show T is a contraction on X. In fact, if u\,U2 £ X, there is 



|T(m)-T(n 2 ) 



2 /"* \/3 
— P / e _s sin — (i — s)(ui — U2)(u\ + u\U2 + it§) (is 
V3 Jo 2 



3 

< 7-75 SUp |ui(t)-ti2(t)| 

2V3 te[o,oo) 

< ^-d{u u u 2 ), 
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which implies 

d(T(m),T(u 2 )) < -^-d(ui,u 2 ). 

So T is a contraction. By Banach fixed point theorem, this theorem is 
proved. □ 

Theorem 2 There is a unique solution y* to the problem 

y"-y' + y = y 3 , -co < x < oo, (2.8) 

(P) 2/(0) = 0, y(oo) = 1, (2.9) 
y(x) > 0, < x < 00. (2.10) 
And y* has infinitely many zeros x n (n = 0, 1,2, ...) 

— oo < . . . < x n < . . . < X\ < xo = 0, 

also 

y(-oo) = 0, \y*(x)\ < ~, 

/or x G (— oo, 0]. 

Proof. The Theorem 2 in || and Theorem 1 imply that problem (P) has a 
unique solution y* . By (p.4|), we see that y* has infinitely many zeros, and 
y{— oo) = 0, |y*(x)| = |n*(-x)ef | < |, for x £ (— oo,0]. □ 

Theorem 3 For any integer n > 0, i/iere is a unique solution y^ to the 
following problem 

y" -y' + y = y 3 ,o < x < oo, 

y(0) = 0,y(oo) = l, 

y has precisely n zeros in (0, oo). 

Proof. By Theorem 2 , it is easy to check that 

y (n \x) =y*( x + x n ) 

is a solution to this problem, where x n is given in Theorem 2. 

Now suppose there is another solution to this problem. Let = xq < 
xi < X2 < ■ ■ ■ < x n be the zeros of y^ in [0, oo). Then y(s) = y( n \s + x n ), 
satisfies (^8|),(^9|),( |2.10| ) with x = s. By Theorem 2 , we have y(s) = 
y*(s), — oo < s < oo, or y( n \x) = y*(x — x n ). Then Xk — x n (k = 0, 1, n) 
are the zeros of y*, y*(xk — x n ) = y^ix^) = 0, satisfying — x n = xq — x n < 
x\ — x n < ■ ■ ■ < Xk — x n < ■ ■ ■ < x n — x n = 0. Because 0, X}.{k = 1, . . . , n) are 
the consecutive zeros of y^ , Xk — x n {k = 0, 1, n) are the consecutive zeros 
of y*. So we see that x n = —x n . Thus y^ n \x) = y*(x — x n ) = y*(x + x n ) = 
y( n \x). So the theorem is proved. □ 
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3 Numerical Results 



In this section, we present numerical computation results for some important 
values related to y*. First of all, as seen in the previous sections, the most 
important number is a*, which defines y* . In the next section, we investigate 
the asymptotics of y*(x) as x approaches to — oo and as x approaches to oo. 
The following values will be used to represent the asymptotics, 



b* = [ e-f CO s^ S (y*( S )) 3 &, (3.1) 

J —OO 4 

c* = f° e-Um—s{y*(s)fds, (3.2) 

J —oo 2 

d* = J™ es^l-y*( s )f- l -(l-y*( s )f}ds. (3.3) 

(1) a* 

In Theorem 2 of |^], we obtained that 

<°*<-I>. 

\ n=l / n=l 

for any iVi > 2, and ^2 > 1- If w e choose large numbers iVi = 170, A^2 = 170, 
the Mathematica shows the following result 

0.16871221576 • • • < a* < 0.16871221594 • • • . (3.4) 

(2) y*(x) and P(z) 

In last section, we have proved that y* (x) is the unique solution to prob- 
lem (P). And it was also shown that y*{x) is the solution of the following 
initial value problem 

y"-y' + y = y 3 , -oo < x < oo, (3.5) 

y(0) = 0, y'(Q) = a*. (3.6) 

So we can approach the function y*(x) by numerically solving the initial 
value problem with properly approximate value of a* given above. The 
Mathematica gives the graph of y*(x). See Figure 1. 

In section 4 and 5, we discussed the function P(z) which is the solution 
to equation PP' - P = z(z - l)(z - 2), z € [0, 1] with P{z) ~ -z as z -» 0. 
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Figure 1: Graph of y*(x)( solution of (P) ). The value of a* = y*'(0) is 
chosen 0.168712 



As mentioned in Q, P(z) gives the relation between y*'(x) and y*(x) for 
x > 0, which is 

y*'(x) = -P(l-y*(x)). 

The graph of P(z) is computed also by Mathematica, and given in Figure 2. 
The numerical computation for P(z) shows that P(l) ~ —0.1687- • •, which 
matches the result (3.4). 
(3) b% c*, and d* 

The values of b*, c*, d* will be used to represent the asymptotics of 
y*(x) as x — > ±oo, and the asymptotics of f(r) as r — > and r — ► oo 
in the next section. We do not have further analytic properties for these 
three numbers since we have not found the analytic representation of y* . 
Computer shows the approximate values of these three numbers based on 
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Figure 2: Graph of P(z) as solution to the equation PP' — P = z(z — 1) [z — 
2),z £ [0, 1] with P(zq) = -z , where z = 10" 40 . 

the numerical computation for y*(x), 

/o fo 
e~i cos Z—3 (y*(s)) 3 ds « -0.0005497 • • • , (3.7) 
-30 2 

f \/3 

c* « / e"t sin -Vs(y*(s)) 3 ds « 0.001939- ■-, (3.8) 
d* « | q 15 ((1 - y*( S )) 2 - 1(1 - y*( S )) 3 ) da « 4.1728 • • • . (3.9) 

4 Asymptotics and Connection Formulas 

We have seen in Sect. 2 that the problem (P) has a unique solution y*(x). 
Now, Let us consider the asymptotics of y*(x) asu -co, and x — ► oo. 
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Theorem 4 The function y*(x) has the following asymptotics 
(i) As x — > — oo, 

y*(x) = Ae% sin {^-x + (pj + O(e^), (4.1) 
with the connection formula 

A = -i\/( a * - b*) 2 + c* 2 » 0.196, (4.2) 

= arctan - - sa 0.0115 sa 0.003757T. (4.3) 

a* - o* 

(m,) v4s x — > 0, there is 

y*(x) ~ a*x. (4.4) 

(mj As x — » oo, 

= l-Be- x + 0(e- 2x ), (4.5) 

luzf/i i/ie connection formula 

B = + d* w 4.90. (4.6) 

Proof. Let us first consider (|4.1| ) . Recall that if we let 

t = — x, u*(t) = e~zy*(x), 
then u* satisfies for t > 

o 

it H — u = e u (t), 
4 

and 

u*(t) = =a*sin^^iH — = / e~ s sin -^-(i — s)u* 3 (s) 

V3 2 V3Jo 2 

= - A a * sin + -| | o °° e" s sin ^(t - s)u* 3 (s) ds + Ofe"*), 

as t —* oo. Direct calculation shows that 

*/ \ 2 ^ . v3 2 ^ y/2> . 3x 

y (x) = —/={a — b )e2 sm —r-x H — =c e^ cos -r-a; + C(e 2 j 
\/3 2 y3 2 

= vle^ sin f-^c + <f> J + O(e^), 



8 



where A and <f> are given by (p~2j),(p~3[). 

The asymptotics ( t4,4| ) is obvious. Finally let us prove ( [4.5|) . Let 



y *(x) = 1 - z(x). 
It's easy to check that z(x) satisfies the differential equation 

z" - z' -2z = -3z 2 + z 3 , (4.7) 

and integral equation 

/ \ 1 — a* o x 2 + a* 

z(x) = e zx H e x + 

3 3 

— \ e- 2s (-3z 2 (s) + z 3 (s))ds-^— e s (-3z 2 (s) + z 3 (s))ds, 

3 Jo o JO 

since z(0) = l,z'(0) = —a*. Using equation ( [L7| ) and Lemma 3.1 (iii) in 
we have 

oo 

e- 2s (-3z 2 (s) + z 3 (s))ds 

oo 

e- 2s {z"{s) -z\s) -2z(s))ds 
z'(0) - z(0) + ((-2) 2 + (-2) - 2) / e- 2s z(s) ds 



Thus 



a — 1. 



= - e -4- [ X e s (-3z 2 (s) + z 3 (s)) ds + 0(e- 2x ) 

3 3 Jo 



2 + a* 



+ d*^j e~ x + 0{e~ 2x ) 



3 

So holds. □ 

Theorem 5 Any solution y(x) to the problem 

y" -y' + y = y 3 , -oo < x < oo, (4.8) 

y(x) — > 0, as x — > —oo (4-9) 

y(x) — ► 1, as x — > oo (4-10) 

has the representation 

y(x) = y*(x-r), (4.11) 

where r is the largest zero point of y(x). Hence the solution to 
(4-&),( 4-~l(\ ) is unique up to a transition of x. 
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Proof. Let y% = y, yi = y', and change equation 



into the system 



y f 2 



2/2, 



V2-yi + yi- 

It is not hard to see that the point (0, 0) in the phase plane is a stable focus 
of the system (see If y(x) is a solution to (fO|) and ( pTIOD , then 

y{x) has infinitely many zeros. Because y(oo) = 1, the zeros of y{x) have 
upper bound, and the largest zero r exists (finite). Let 

2/1 0*0 = y(x + r). 



Then yi(x) solves (P). By Theorem 2 , y\(x) = y*(x), and then (4.11) holds. 



□ 



The solution to problem (P) is clear now. Let us come back to the 
original equation (1.1). 

Theorem 6 Any solution to \1. f\) , and has the representation 



f(r) = y*(log 



ro 



(4.12) 



where ro is the largest zero of f(r). With a*, b*, c* and d* defined by 



(1.4), ^J\ ), (3.i) and {3S} respectively, f(r) has the following asymptotics. 
(i) As r -> 0, 



S(r) = A ( ^) 2 sin log 1- + ^ + 0(A), 



where again 



A 



ro 
2 

71 



\J(a* - b*) 2 + c* 



0.196, 



arctan 



a* — b* 



0.0115 « 0.00375vr. 



(ii) As r — > ro, 
(mj As r — > oo, 



f(r) ~ a* log—, 



where 



B = — — + d* « 4.90. 



(4.13) 

(4.14) 
(4.15) 

(4.16) 

(4.17) 
(4.18) 



Proof. By Theorem 4 and Theorem 5 , using r = e^ro = e r ,log^ 
x — t, /(r) = y(x), direct calculations show the results. □ 
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